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This  final  report  was  prepared  by  the  University  of  New  Mexico, 
Albuquerque,  New  Mexico,  under  Contract  F29601-84-K-0045,  Job  Order  ILIR8423. 
David  H.  Hilland  (NTA)  was  the  Laboratory  Project  Officer- in-Charge. 

When  Government  drawings,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  in  connection  with  a  definitely  Government-related 
procurement,  the  United  States  Government  incurs  no  responsibility  or  any 
obligation  whatsoever.  The  fact  that  the  Government  may  have  formulated  or 
in  any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not 
to  be  regarded  by  implication,  or  otherwise  in  any  manner  construed,  as 
licensing  the  holder,  or  any  other  person  or  corporation;  or  conveying  any 
rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention 
that  may  in  any  way  be  related  thereto. 

This  report  has  been  authored  by  a  contractor  of  the  United  States 
Government.  Accordingly,  the  United  States  Government  retains  a  nonexclu¬ 
sive,  royalty- free  license  to  publish  or  reproduce  the  material  contained 
herein,  or  allow  others  to  do  so,  for  the  United  States  Government 
purposes . 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  and  is 
releasable  to  the  National  Technical  Information  Services  (NTIS).  At  NTIS, 
it  will  be  available  to  the  general  public,  including  foreign  nations. 

If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing 
list,  or  if  your  organization  no  longer  employs  the  addressee,  please  notify 
AFWL/NTA,  Kirtland  AFB,  NM  87117-6008  to  help  us  maintain  a  current  mailing 
list. 

This  technical  report  has  been  reviewed  and  is  approved  for 
publication. 
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We  looked  at  32  fly-by  (Ref.  1)  data  sets  from  3  different  locations 
(fuselage,  left  wing,  and  right  wing)  for  4  different  directions  of 
incidence  (nose  on,  left  wing  on,  tail  on,  and  right  wing  on).  The  sampling 
rate  was  10  ns  for  every  data  set. 

Out  of  32  data  sets,  only  6  data  sets  from  the  fuselage  show  consist¬ 
ency  in  the  first  3  pole-pairs  (Fig.  l).  Plots  of  these  six  data  sets  are 

shown  in  Fig.  2.  Plots  of  the  Fourier  transform  magnitude  of  each  waveform 
are  shown  in  Fig.  I.  Note  that  there  is  a  DC-bias  in  some  waveforms.  This 
was  corrected  before  computing  a  spectrum.  We  also  omitted  the  first  few 
time  samples  in  an  attempt  to  remove  some  of  the  incident  field. 

Before  applying  the  pole-finding  algorithm,  each  waveform  was  passed 
through  a  second  order  digital  Butterworth  lowpass  filter  with  a  20  MHz 
cutoff  frequency.  The  number  of  data  points  used  was  90.  The  approximate 
iterative  preprocessing  algorithm  (AIPA)  (Ref.  2)  for  multiple  measurements 
was  applied  for  six  data  sets.  The  first  couple  of  pole-pairs  were 
successfully  extracted  by  processing  all  six  data  sets  at  the  same  time. 

Pole  #1:  =  -6.715  x  10s  ^  =  6.803  Wz 

Pole  #2:  a2  =  -8.602  x  10s  f2  =  11.509  Wz 

The  number  of  iterations  required  for  convergence  was  less  than  10. 

The  third  pole-pair  was  extracted  by  processing  data  sets  for  left  wing 
on  and  right  wing  on.  When  the  order  was  6,  the  result  was  as  follows: 


1.  Lee,  K.  S.  H.  et  al.,  “Preliminary  Comparison  of  F-106B  Responses  to 
Natural  Lightning  and  Simulated  NEMP,"  Air  Force  Weapons  Laboratory, 
Heart  Conference,  Summer  1985. 

2.  Park,  S.,  and  Cordaro,  J.  T.,  Preprocessing  Techniques  for  the 
Estimation  of  Parameters  from  Noisy  Transient  Data,  AFWL-TR-85-102,  to 
be  published. 
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d1  =  -8.263  x  10s 


a2  =  -7.334  x  10s 


a3  =  -6.423  x  10s 


fj  =  6.043  MHz 


f2  =  11.625  MHz 


f,  =  15.457  MHz 


When  the  order  was  8,  the  result  was  as  follows: 


-7.653  x  10s 


-8.082  x  10s 


-6.656  x  10s 


-3.123  x  10s 


=  6.091  MHz 


11.567  MHz 


=  15.281  MHz 


=  19.914  MHz 


Because  the  fourth  pole-pair  was  not  numerically  stable,  we  selected  the 
third  pole  as: 


Pole  #3:  a3  =  -6.656  x  10 


fs  =  15.281  MHz 


Using  the  algorithm  presented  in  Chapter  7  (Ref.  2),  the  amplitude  of 
the  normalized  coupling  coefficients  for  the  first  three  poles  were 
calculated.  These  are  shown  in  Fig.  3. 

Theoretically,  the  amplitude  of  the  normalized  coupling  coefficients 
for  90  deg  (left  wing  on)  and  270  deg  (right  wing  on)  should  be  the  same 
(by  symmetry);  but  the  plots  in  Fig.  3  do  not  show  exact  symmetry — although 
plot  (c)  is  close.  Only  one  data  set  for  each  wing  on  orientation  was 
available.  Also,  the  transform  amplitudes  for  poles  1  and  2  at  90  deg  and 
270  deg  are  quite  different  as  can  be  seen  in  Fig.  1(c)  and  (f).  As  a 
result,  the  expected  symmetry  was  not  found.  This  could  be  due  to  experi¬ 
mental  error:  for  example,  the  aircraft  may  not  have  been  positioned 
exactly  wing  on  for  these  particular  shots. 
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As  shown  in  Chapter  4  (Ref.  2),  the  shorter  the  sampling  rate,  the 
smaller  the  percentage  error  in  the  pole  estimate.  Data  sets  with  sampling 
rate  of  5  ns  should  give  better  results. 


The  iterative  preprocessing  algorithm  (IPA)  (Ref.  2)  was  also  used  to 
find  poles,  but  the  results  were  almost  identical. 


